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Abstract 
> 

We use a WKB approximation to establish a relation between the wave- 
front velocity in a strongly coupled theory and the local speed of light 



in a holographic dual, with our main focus put on systems with Lifshitz 
scaling with dynamical exponent z. We then use Einstein equations to 
relate the behavior of the local speed of light in the bulk with the null 
energy condition (NEC) for bulk matter, and we show that it is violated 
for Lifshitz backgrounds with z < 1. We study signal propagation in 
S^ the gravity dual and show that violations of the NEC are incompatible 

with causality in the strongly coupled theory, ruling out as holographic 
models Lifshitz backgrounds with z < 1. We argue that causality vi- 
olations in z < 1 theories will show up in correlators as superluminal 
modes and confirm this for a particular example with z = 1/2. Finally, 
as an application, we use z < 1 solutions to uncover regions of the 
parameter space of curvature squared corrections to gravity where the 
NEC can be violated. 



1 Introduction 

The AdS/CFT correspondence fil has been extensively used as a tool to extract properties 
of strongly coupled systems. In its usual formulation the strongly coupled dynamics of a 
large- iV gauge theory is extracted from a gravitational theory in a higher dimensional space. 
Recent developments are directed towards the development of similar techniques for strongly 
coupled critical points appearing in condensed matter systems, a nice introduction to the 
subject can be found in refs. (2l. Although critical points show some kind of scale invariance, 
Lorentz symmetry is usually broken so the time coordinate can scale differently to the space 
coordinates 

t^XH, x^Ax. (1.1) 

Here z is known as the dynamical critical exponent. In spatially anisotropic systems where 
rotational invariance is broken one or more spatial coordinates can also scale differently to 



the rest. Systems with dynamical scaling ( 1.1 ) have been studied for a long time in condensed 
matter theory J3J. 

In holographic models with dynamical scaling the starting point is to construct a ge- 
ometry whose isometries and local symmetries map to the global symmetries of the critical 
point of interest, most importantly to the scaling properties. Examples where this was first 
done are refs. [4|j6]. So far theories with known holographic duals are still deformations of 
large- N gauge theories, but for properties that do not depend strongly on the microscopic 
details of the theory one wants to study, these holographic constructions could produce good 
qualitative results that would be very difficult to obtain using other methods. 

In a relativistic theory the dispersion relation for massless particles is such that the 
frequency is proportional to the momentum u = ck. Theories with broken Lorentz invariance 



z 7^ 1 (1.1) can have different dispersion relations at small frequencies, so in general an 
ultraviolet completion would be needed in order to make them compatible with causality. 
Consider, for instance, a mean field description of massless scalar fluctuations around a 
critical point with dynamical exponent z. The effective Lagrangian would be 

£=(9 t 0) 2 -c 2 ^- 1 V(-^r0, (1.2) 

where £ has units of length and c is the speed of light. The dispersion relation reads 

u? = j 2 {ikf z - (1-3) 

It follows that the phase velocity is 

v ph = ^ = c(£k) z - 1 . (1.4) 

We could think of l/£ as a cutoff in k, at k = l/£ fluctuations reach the speed of light and 
for larger values the theory should have a relativistic description. Although strictly speaking 
the phase velocity does not have to be smaller than the speed of light, in the limit u — > oo 
(so k —?• oo) the phase velocity becomes the wavefront velocity, that should not be larger 
than the speed of light c. The limit where relativistic effects can be neglected corresponds 
to c — > odj and £ — >• while keeping k = c£ z ~ x fixed. 

1 Strictly speaking we rescale the time direction so the speed of light seems to grow. 



In the discussion above we have made the implicit assumption that z is an integer, 
otherwise the term with spatial derivatives would not be local. Fractional values of z may 
appear in more complicated situations. We should also remark that if z < 1 the previous 
argument fails, the phase velocity diverges when k — > and does not when k — > oo. Clearly 
a theory with dynamical exponent z < 1 cannot be an infrared fixed point, although it is 
not ruled out that such a theory describes an intermediate range of scales in some system, 
as for instance this could be the case for fully developed turbulence [7]. 

The metric proposed for a holographic description of a critical point with dynamical 
exponent z is |5l[8] 

ds2 = ^{-^=r ) + dr2 + d ^)- (1-5) 

Here L is the radius of curvature of the metric and there is a conformal boundary at r = 0. 
One can apply a holographic recipe adapted to these spaces to compute two-point functions 
in the dual theory (5l[9]. The results for scalar operators agree with expectations from scale 
invariance, in particular they depend on the combination u/(nk z ). 

Now let us make the following observation: the local speed of light at a fixed value of the 



radial coordinate has a dependence on 1/r that is the same as the phase velocity in (1.4) 
with k, namely 

c {r)=c^ z - l) r- {z - 1) . (1.6) 

This observation fits well with our intuition of how holography should work, the radial 
coordinate is associated with different scales, and radial slices to the field theory at those 
scales. We will make this argument more rigorous and actually show that the local speed of 
light is always related to the phase velocity at small wavelengths even in situations where 
there is no exact scale invariance. 

The next point we will address is what is the difference between z > 1 and z < 1 from 
the gravitational perspective. Although both cases have different singular behavior, it is not 
a priori clear why one choice would be 'better' than the other. We will show that z < 1 
backgrounds are incompatible with causality in the holographic dual, and that this will be 
true in general for any geometry with a local speed of light decreasing towards the boundary. 
We will show that violations of causality in this sense are produced by matter that violates 
the null energy condition (NEC). 

We also study two-point functions of scalar operators computed following the holographic 
recipe. We write the equation of motion for scalar fields as a Schroedinger equation to argue 
that the qualitative behavior in z < 1 and z > 1 geometries is completely different. For z < 1 
and any finite value of k, the potential is confining so the holographic two-point function 
should show a discrete set of poles. When k — > the poles merge forming a branch cut. We 
give an explicit expression for z = 1/2 that confirms our arguments. 

The current paper is organized as follows. In Sec. [2] we study the relation between local 
speed of light and phase velocity in the limit where the phase velocity becomes equal to the 
wavefront velocity. We then argue about causality in Lifshitz backgrounds for z > 1 and 
z < 1 and find very different structure in those cases. Sec. [3] is devoted to the calculation 
of two-point functions for scalar operators using holographic techniques. We find poles in 
the propagators for different critical exponents and analyze spectra of scalar perturbations. 
Having established the importance of the NEC for the holography we can use it to constrain 



some known theories. In Sec. [4] we study higher derivative gravity theories which are also 
used in holographic constructions and determine a region where Lifshitz solutions with z < 1 
exist and hence violations of the NEC are possible. In Sec. [5] we present our conclusions and 
discuss some open questions. 

2 Null energy condition and causality 

We will use a simple setup with a scalar operator in a d + 1 dimensional field theory with a 
D = d + 2 gravity dual. The holographic description of the vacuum is given by a background 
metric and fields in the gravity dual. Assuming rotational invariance in the spatial directions, 
we can write the metric as 

ds * = du 2 + e 2A(u)(_ e 2B(u) dt 2 + ^2) ^ ^ 

where t, x correspond to the time and space coordinates in the field theory and u is the 
holographic radial coordinate. In this coordinate system the boundary is located at u — > oo. 

We can now introduce a source for the operator for a finite interval of time and let the 
system relax. If the perturbation is small enough, so linear response theory holds, the final 
state will simply consist of some scalar modes propagating through the vacuum. One can 
then expand in plane waves to study the dispersion relation of the scalar modes. The phase 
velocity is given by the ratio of frequency and momentum v p h = ui/k. In the u — > oo limit, 
the phase velocity becomes the wavefront velocity v wi , that for a relativistic theory should 
be smaller or equal to the speed of light v w f < 1. In a non- relativistic theory the velocity 
can take any value, but it should remain finite at large frequencies if the theory is to have a 
relativistic completion, as has been discussed in the introduction. 

In the holographic description the states created by a scalar operator correspond to 
classical normalizable solutions of a dual scalar field. We have just to consider the quadratic 
part of the action 

S = - I d d+2 x^ (d M ® d M $ + m 2 $ 2 ) , (2.2) 

where the mass m depends on the scaling dimensions of the operator. Using a plane wave 
ansatz $(£, x, u) = e~ luJt+lkx (j)(u) the equations of motion read 

0" + ((d + I) A! + B')4>' + e- 2A - 2B uj 2 (j) - e- 2A k 2 $ - m 2 = . (2.3) 

It is useful to rewrite this equation as a Schrodinger equation. First we define a new coor- 
dinate 

p '(u) = _ e -A{u)-B(u) _ ^^ 



Then we can rewrite (2.3) as follows (derivatives with respect to p will be denoted by dots) 

4> + dA<j) + u 2 <P - k 2 e 2B <P - m 2 e 2A+2B <P = . (2.5) 



Now define <fr — cip with 
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The equation (2.5) then becomes 



- 4> + V{p)ip = u 2 ^ , (2.7) 

which is a Schrodinger equation with u 2 playing the role of the energy. The potential iqj 



Via) = k 2 e 2B + m 2 e 2A+2B + -A 2 + -A . (2.8) 

\r, 4 2 

In order to make the construction more explicit let us consider the Lifshitz metric. Then in 



the notations of (2.1) we have 



A(u) 



u 

L 



B{u) = {z-l) 



(2.9) 



and the variable p = Le zu ^ L /z, such that p — > is the boundary limit and p — > oo is the 



horizon limit, this is equivalent to doing the change of coordinates p = r z /z in (1.5). We 
rescale the time coordinate so u 2 —¥ u 2 /n 2 . Then the potential (2.8) reads 



T _. . d 2 + 2dz + Am 2 L 2 ,o/pz\- 



-2 



(2.10) 



If we use the relation between the mass and the scaling dimension of the dual operator 
m 2 L 2 = A(A — d — z), the potential becomes 



V( P ) 



A 



d + z 



^2 



-{zp? /z 



(2.11) 



The potentials (2.10) for d = 2 and m 2 = are plotted in Fig. n] for several values of z 



The first term in (2.10) is singular near the boundary, while the behavior of the second term 



depends on the value of z. By simple inspection of the potential we can see that solutions 
have a different qualitative behavior depending on whether z > 1 or z < 1. If z > 1, the 
potential decays as p — > oo, so the solutions become plane waves ip ~ e ±tU)p and the spectrum 
of normalizable solutions is continuous. For z < 1 and k ^ the potential has a barrier at 
large values of p. There is a minimum at 



2r2^ 



L ((l-z)(d 2 + 2dz + 4m 2 L 
4z{kL) 2 



(2.12) 



Notice that for z = 1 and k 2 > u 2 , solutions are either exponentially growing or decreasing. 
Only the latter can belong to the physical spectrum of the theory, that should be normalizable 
(or delta-normalizable) over the entire range of the p coordinate. For smaller values of z 
(z < 1) and k y^ the barrier is steeper, so normalizable solutions will be suppressed as well. 
At p = some negative values of the mass allow two possible normalizable solutions (when 
the coefficient of the 1/p 2 term is between 3/4 and —1/4), but in holographic applications 
each choice corresponds to a different boundary theory [12] . Imposing normalizability at 
both ends leads to a discrete spectrum. 



2 Similar analysis was performed in refs. 
bulk were considered. 
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where braneworlds with broken Lorcntz invariance in the 



V(p) 




Figure 1: Quantum mechanical potentials for z = 0.2,0.4,0.6,0.8 as functions of the radial coordi- 
nate p. The plots are ordered by their slopes at large p, the steepest curve corresponds to z = 0.8. 
Spatial momentum k is set to unity together with the bulk curvature scale L. 



Phase velocity from the WKB approximation. Consider the limit of large momentum 
k — > oo and large frequency u —¥ oo. The term proportional to k 2 in the potential (2.8) 
dominates, except very close to the boundary p = 0. In this limit we can use the WKB 
approximation to the solutions. Notice that e B is the function that gives the local speed of 
light at a fixed value of the radial coordinate. If we assume that it decreases, then for large 
values of p the potential will be smaller than the energy u 2 > V(p) and the WKB solution 
will be oscillatory 



V>(p) 



C\ 



n(u 2 -V(p)) 



1/4 



ifdp^w*-V(p) 



+ 



(■■2 



n(uj 2 -V(p)) 



1/4' 



ifdpyJuP-VUi) 



(2.13) 



However, for values of p very close to the boundary V(p) ~ m 2 e 2A > u 2 . This means that 
there is a turning point po where V(po) = u 2 . In the interval where V(p) > u 2 and the k 2 
term in the potential still dominates, the WKB solution will be an exponential 

^1 CdnJv(n\-ufl , d 2 



4>{p) 



J'dpy/Vifi-J* 



+ 



f p dp^V(p)-^ 



(2.14) 



/Jr (V(p) - u 2 ) 1/4 ^{V{p)-u 2 ) 1/4 

As p — > the potential is too steep V(p) ~ 7/p 2 for the WKB approximation to be valid (see 
for instance [13] for a treatment of both the WKB approximation and the 1/ p 2 potential). As 
long as 7 > —1/4 (as we will assume it is the case), one can simply do a Frobenius expansion 
at p = 0, pick the normalizable solution and then match to the WKB solution (2.14). The 



solution at large values of p (2.13) can be fixed by matching both WKB solutions at the 



turning point p with Airy functions. 

We are not really interested in finding the solution, but rather in the value of the turning 
point in the limit u — > 00 we are taking. The condition V(po) = u> 2 gives us 

to 



Vwf - V ph 



Mpo) 



k 



(2.15) 



We see that in this limit the wavefront velocity becomes equal to e B ( p0 \ which is the local 
speed of light at the turning point. Therefore plane wave states created by a scalar operator 
in the field theory have wavefront velocities that are equal to the local speed of light in the 
holographic dual. This confirms the observation we made in the introduction in relation to 
the local speed of light in Lifshitz geometries. 

Growing versus decreasing speed of light. In the previous paragraph we assumed 
that the local speed of light on a radial slice decreases when the radial position is moved 
towards larger distances from the boundary, and we have assumed m 2 > as well. In 
principle this argument can be generalized to cover more cases, the potential will be slightly 
more complicated and whether the solution is oscillating or an exponential on different radial 



intervals can vary. However, the matching condition (2.15) that gives the relation between 
phase velocity and local speed of light would be the same in the u> — > oo, k — > oo limit. 

It is interesting to study the difference between geometries with a local speed of light that 
grows towards the boundary and geometries where it decreases. The Lifshitz metric gives 
examples for both, depending on whether z > 1 (growing) or z < 1 (decreasing). We have 
already commented in the introduction that field theories with dynamical exponent z > 1 or 
z < 1 have qualitatively different dispersion relations, in the former the phase velocity grows 
with momentum, while in the latter it decreases. We have just shown how this is reflected 
in the local speed of light in the geometry. 

A consequence of this difference in dispersion relations is that the boundary structure of 
z > 1 and z < 1 geometries is quite different. We shall recall here the Lifshitz metric 

ds>=^(-^+dr* + dA. (2.16) 



nr> A \ /yAyZ L ) 

For z > 1 we do the change of variables R = r z , t — > t/z, x — > x./z 

ds 2 = -^- (-n 2 dt 2 + dR 2 + R 2 ^ z d^ 2 ) . (2.17) 

z 2 R 2 

The conformal boundary is at R — > 0. Notice that for z > 1 it is timelike and one- 
dimensional, or in other words it is along the time direction. For a general value of z it 
is also singular. For instance, in the particular case of z = 2 there is a conical singularity. 
When z — > oo the geometry is AdS2 x M. d , and the boundary is regular. 
For z < 1, we can write the metric in the original coordinates as 

ds 2 = ^ (dr 2 + dx 2 - r 2 - 2z K 2 dt 2 ) . (2.18) 

In this case the conformal boundary r — y is a <i-dimensional surface. In general the 
boundary is also singular, for instance when z — 1/2 the singularity is conical. For z — 1, 
when the geometry is AdSd+2 the boundary is regular and it also includes time. 

Notice also that for z > 1 the slope of null geodesies in the (t, r) plane is such that they 
are orthogonal to the boundary, while for z < 1 they are tangent, so the singularity is null 
instead of timelike. Indeed, for z < 1 the slope becomes infinite at the boundary but it is 



reached in a finite time (for z > 0), the equation for a null geodesic is 

dt r z ~ l 

* = — iM = - (2 - 19) 

which leads to 

t(r) = ^ . (2.20) 

ZK 

The presence of a singularity at the boundary, or equivalently at large scales in the dual 
theory, suggests that an ultraviolet completion is necessary^} A possibility is to cut the 
geometry at a small value of r, r = P. <C L. For Lifshitz backgrounds with z > 1 we could 
identify £ with the cutoff in the non-relativistic theory and KJl z ~ l with the speed of light 



of the relativistic completion. In some cases 15 the geometry is only Lifshitz close to the 



horizon and becomes an AdS space close to the boundary, so the ultraviolet completion 
has already been included in the description. A different approach is to try to define the 
holographic theory without changing the asymptotic boundary behavior or introducing a 
cutoff, as in ref. |16| where the notion of conformal boundary is extended to anisotropic Weyl 
transformations. A well defined prescription requires to constrain the allowed fluctuations 
in the bulk, in some cases this implies constraining the sources of the dual theory [l7|. A 
formal treatment clarifying these issues is certainly desirable. For our purpose it will be 
enough to introduce a cutoff, since the results we will obtain are independent on how the 
ultraviolet theory is defined. 

A natural question is whether a consistent holographic description requires further condi- 
tions. For instance, we have commented in the introduction that Lifshitz theories with z < 1 
would have a bad infrared behavior. However, the boundary analysis of the holographic 
dual does not seem to be so helpful here. Both z > 1 and z < 1 geometries have a singular 
boundary, and we would introduce a cutoff to get rid of it. In the following we will use 
causality arguments to show that even with a boundary cutoff the z < 1 geometries are not 
good holographic duals. 

Causality from shock waves. We will now perform an analysis in the spirit of the 



conformal colliders thought experiments of Hofman and Maldacena 18 . We will introduce 
a source in the field theory localized in time and in one of the spatial directions x. This 
will produce a planar perturbation propagating along x. In the holographic description the 
source is a boundary condition that will produce a shock wave-like perturbation propagating 
along null geodesies in both the radial and the spatial direction. This corresponds to both 
the propagation of the wave and its spreading to larger wavelengths. The shock wave will be 
a source of radiation of gravitational fields that will then propagate along the radial direction 
to the boundary, producing a front of radiation that can be interpreted as the front of the 
perturbation in the dual theory. A similar picture will arise if one uses a probe dragging 



string, as has been done to study jet physics in holographic duals 19 . An alternative way of 
measuring the effect of the shock wave is to have static probes consisting of strings extended 
in the radial direction. When the shock wave crosses a probe, it will produce a perturbation 



3 This does not mean that the bulk theory is not well defined in the presence of singularities. Although 
classically there is a singularity, the motion of quantum test particles could be well behaved 114] . 



that will propagate towards the boundary along radial null geodesies producing a signal that 
can be interpreted as the position of the front. We have illustrated this picture in Fig. [2} 



boundary 




shock wave 



Figure 2: A source at the boundary produces a shock wave propagating through the bulk. When 
the shock wave crosses a static probe a pulse of radiation is sent to the boundary. In the holographic 
dual the front of radiation produced by the source is determined by the the time and position of 
the emitted bulk radiation when it reaches the boundary. 



We can find the null geodesic^j by solving the variational problem with Lagrangian 



£ 



K 2 i 2 x 2 + r 2 

+ „ — 



.-2z 



(2.21) 



and imposing the constraint £ = 0. Since £ does not depend on t or x, we can introduce 
two integration constants that solve t and x in terms of r 



t 



Er 



2- 



x = Pr* 



(2.22) 



If we now solve the constraint 






2 /J 5 y.(,-i) p2 



/■v 



Then the geodesic equations in terms of r read 



dt 
dr 



Er 2(z-1) 



2 / £2 r 2(z-l) _ „ 2 



dx 
dr 



E 2 T 2(z-l) 



-P 2 



(2.23) 



(2.24) 



Let us assume that the initial point for the geodesic describing the shock wave is at r = £, 
t = 0, x = 0. In order to have a sensible solution we need the momentum to be timelike 

p 2 = 2 P 2 > 0. We see now a clear different qualitative behavior between z > 1 and 

z<l. 

4 Recall that in 11 null geodesies were used for qualitative studying localization of field fluctuations in 
the bulk. 



For z > 1 the argument inside the square root in (2.24) is always positive, so the geodesic 



extends to all values of r. The asymptotic behavior at large r is 

r z kP 



t 



~ > oo . 

ZK 



X 



(2 - z)E 



r 2 z + x 



(2.25) 



For the special value z = 2, x grows logarithmically with r. We see that for z > 2 signals at 
the boundary can reach only a finite distance that depends on the ratio P/E < £^ z ~^ / k. 



For z < 1 the argument inside the square root in (2.24) becomes negative at a finite value 

c£ z ~ x and assuming P > 0, 



of the radial coordinate, when r 



1-2 



E/(k\P\), or using k 



To 



E 
cP' 



(2.26) 



Although the slope diverges, the values of both t and x approach a finite value, so the null 
geodesic is tangent to the slice of constant r at that point, and it will bounce back towards 
the boundary. In the above formula tq represents the right turning point of the quantum 



mechanical potential (2.11). 



The local speed of light at r = £ is c, that is also the speed of light of the dual theory, we 



can compare this value with the average velocity of the shock wave. Using (2.24) and (2.26) 



the space interval that the shock wave has traveled when it comes back to the boundary is 



Ax 



cP 



'■o 



dr 



(r/£) 



1-2 



E 2 



(?)' 



■''o 



dr 



(r/r ) 



1-2 



1- 



2(1-2) 



(2.27) 



Similarly, the time interval is 

To 



cAt 



dr 



(r/£) 



2-1 



1- 



C 2 P 2 

E 2 



(5) 



2(1-2) 



cP 



■'•() 



dr 



(r/r ) 



2-1 



1 - 



2(1-2) 



(2.28) 



Let us now assume that tq 3> £, so the turning point is located far apart from the boundary. 
Then we can take the limit £ — > in the above integrals, and after changing the variables 
P = r / r o w e get 



Ax 



2r / dp 

'0 



P 



1-2 



Trr 



p 



2(1-2) 



cAt 



cP 

2 B r ° 



„2-l 



dp 



P 



2(1-2) 



2r n 



cP 

2 B r ° 



2-2 
2(1-2) 



r 



1 



2(1-2) 



2(1-2) 



(2z-l)T 



22-1 
2(1-2) 



Then the shock wave velocity is 



c 



Ax 
cAt 



E 
'^P 



(2.29) 



(2.30) 



(2.31) 



Since r 3> £, due to (2.26) the latter ratio is bigger than unity, thus the shock wave travels 



faster than light signals at the boundary v$ > c. One can show that the front of radiation 
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will coincide with the shock wave at the boundary, so its average velocity is also larger 



than the boundary speed of light. As was argued in ref. 20 for asymptotically AdS spaces, 
whenever there is a null geodesic returning to the boundary there are singularities in the 
correlation functions that correspond to poles in the Fourier transform. Similar arguments 
should apply in this case, but the poles will correspond to superluminal modes, we will 
show this explicitly in Sec. [3j Superluminal propagation of this kind is incompatible with a 
holographic interpretation of a causal theory, so I must be an infrared cutoff. Notice that in 
this case the non-relativistic limit in the field theory will be c — > oo, £ — > oo and k fixed, but 
an ultraviolet cutoff is still needed in order to have a relativistic completion. As we expected 
from field theory arguments, the z < 1 theory can only be well defined in an intermediate 
range of scales. 

The argument for Lifshitz geometries can easily be extended to any geometry with a local 



speed of light that decreases monotonically towards the boundary. Using the metric (2.1), 
the variational probe for null geodesies has the Lagrangian 

C = -e 2A+2B i 2 + e 2A x 2 + u 2 . (2.32) 

We can follow the same steps and introduce the conserved quantities E and P 

i = Ee- 2A - 2B , x = Pe~ 2A 1 (2.33) 

which is followed by 

u = e~ A VE 2 e~ 2B - P 2 . (2.34) 

Then the geodesic equations in terms of the radial coordinate are 

dt Ee- A ~ 2B dx Pe~ A 



du y/E 2 e- 2B - P 2 du ^E 2 e~ 2B - P 2 

The local speed of light is simply e B , so the argument of the square root will become negative 
at a finite value of u for a null geodesic starting at the boundary. The arguments we have 
used for Lifshitz geometries hold for more general cases, geometries with such behavior all 
the way to the horizon are not sensible holographic duals. 

Speed of light and the null energy condition. The null energy condition for Lifshitz 
geometries was studied in ref. |8|. There it was shown that backgrounds with z < 1 violate 
this condition, as opposed to backgrounds with z > 1. A natural question is then if the 
behavior of the local speed of light is related to the null energy condition so a consistent 
holographic description requires it. We will show that this is indeed the case. 



We will use the metric (2.1) for the derivation. The NEC for a field theory implies that 



its energy-momentum tensor is semi-positive definite on the light cone 

T^CC>0, (2.36) 

for an arbitrary null vector £ M . For example for the perfect fluid this condition is transformed 
to p + p > 0, where p and p are pressure and energy density of the fluid, or if one introduces 



11 



the equation of state p = up, the condition implies u > — 1, becoming exactly uj = — 1 for 
the cosmological constant. 

The energy-momentum tensor enters in the Einstein equations 

Rfj.v ~ t;9llvR + Ap Miy = T^v , (2.37) 

where Rn U is the Ricci tensor and we have set the Newton's constant to be equal to one. Using 



the Einstein equations we can recast (2.36) in terms of the metric (2.1) and its derivatives 



only without specifying the matter content of the theory. As it was shown in |8|, for a 



diagonal energy-momentum tensor in presence of the spatial SO(D — 2) isotropy (2.36) is 
equivalent to 

Rl-R x x <0, Rl-R u u <0, (2.38) 

where the components of the Ricci tensor in the D dimensional bulk read 

R\ = -B" - DAB' - B' 2 - A" -(D- l)A' 2 , 

R x x = -A'B' - A" -(D- I) A 12 , 

Rl = -B" - (A' + B'f -{D- I) A" -{D- 2)A' 2 , (2.39) 



where the primes stand for the derivatives with respect to u. The first inequality (2.38) leads 
to 

B" + B'{B' + (D-1)A') > 0. (2.40) 

Let us first see what happens when the null energy condition is satisfied. We define a function 
of the radial coordinate C(u) such that 

B' = Ce-( D -^ A - B . (2.41) 

Then the condition reads 

C , e -(D-1)A-B> Q} (2>42) 

or merely C > 0. The derivative of the local speed of light is 

(e B ) , = B'e B = Ce- (D ~ 1)A J (2.43) 

whose sign is the same as the sign of C. Now we can study the different possibilities. First, 
if C > for some value of u = u*, then C > for all u > u* and the speed of light is 
monotonically increasing with u. If instead C < for some value of u = m*, then there will 
be a value uu > w* such that C{uh) = and we will go back to the first case. There could 
be a fine-tuned situation where uh — > oo, but it is unclear whether that is possible or not. 

If the null energy condition is violated, we can define the function C(u), but now C < 0. 
Then, if C < for some value of u = u* the local speed of light will be monotonically 
decreasing for u > u*. If, on the other hand, C > for some value of u = w*, there will be 
a value uh > «* such that C(uh) = 0. Again, there could be a fine-tuned situation where 
u H —7- oo. 

We have then proven that in generic situations the null energy condition is necessary 
in order to have a consistent holographic description. In particular, for the Lifshitz metric 

12 



(1.5) the NEC is satisfied if z > 1, only these values of z are suitable for a holographic 



construction. Notice that there could be causality issues of the type we have described when 
the space is asymptotically AdS (z = 1) if there is matter in the bulk that violates the 
NEC 
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For our considerations we do not need the second inequality from (2.38 ) since the behavior 



of the local speed of light as a function of u is governed by tt and xx components. In refs. 22 



the second condition from (2.38) was used to model the RG running of central charges in a 



conformal field theory using a holography dual. 



3 Scalar two-point correlation functions. 

We now study how the value of the dynamical exponent affects the correlation functions of 
operators in the field theory. We will consider scalar operators for simplicity, but we expect 
that the main qualitative features will be generic for operators of different spin. A scalar 
operator Oa of scaling dimension A in the field theory is dual to a scalar field <p of mass 
m 2 L 2 = A(A — d — z) in the Lifshitz geometry (1.5). The scalar correlator can be computed 



following the usual procedure of evaluating the action on a classical solution. The action 
is in general infinite, so it needs to be properly renormalized through the introduction of 
boundary counter-terms. We are not interested in the ultra-local behavior of the correlator, 
so we will ignore this issue and just keep the finite piece of the action. For a more complete 
discussion see ref. 23 . Euclidean correlation functions were originally computed for the 



z = 2 case in ref. (5 . 

Given a plane-wave solution 0(t, x, r) 
holographic two-point function is 



e luJt+lkx tpuj,k{r) of the equations of motion, the 



G 2 (w,k) = -lim v /= ^^"^L,k( r )^,k( 



\r=£ 



(3.1) 



The usual 



We have to specify boundary conditions for the solution in the r — ?• oo limit, 
identification is that an ingoing wave corresponds to a retarded Green's function 24 , however 
for z < 1 it is not always possible to choose an ingoing condition. Instead, we will fix 
the solutions to be regular and when an ingoing solution is possible, to be the analytic 
continuation of a regular Euclidean solution. Therefore we will be computing Wightman 
correlators. 

G 2 (w,k)= [ dtd d xe tiVt - lkx (O(t,x)O(0)) . (3.2) 



The equation of motion for a scalar field in (1.5) is 

+ d 



<p"~ 



-(ft + —r 2{z - 1} (p - k 2 (p - 

k 2 ' ' r 



i 

2^ 







(3.3) 



where k = |k|. Notice that k has a mass dimension that depends on z, [k] = 1 — z. 

Let us now compute the two-point functions for three special values, z = 2, z = 1 and 
z = 1/2 l^j For illustration purposes, it is enough if we consider massless fields A = d + z 
and fix the number of spatial dimensions to d = 2. The solutions satisfying the appropriate 



3 The cases z = 0, z — > oo can be found in ref. |9|. 
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boundary conditions are of the form </?(r) = $(r)/$(£), where $(r) is the continuation of 
the regular Euclidean solution. We list our results in the following 



z=2 (Lifshitz) There is a subtlety here, the solution depends on whether we want to 
extend the correlator over the upper half of the frequency complex plane or over the lower 
part. The solution for Imw < is 



*<M = e- >-{/(---— ,-1,— 



(3.4) 



where U(a, b, c) is a confluent hypergeometric function. Meanwhile, the solution for Imw > 

is 



<K(r) = e 2k 



L x 



iujt 



hi 



ie 4^ / 1 ink 2 iur" 

•■'3 ; fcM^l - ^"!^"'" 1 ' 



■"■ I 2 ~ 4<W 

where L^(x) stands for a generalized Laguerre polynomial 



K 



(3.5) 



z=l (AdS) 



$0) 



-VP 



1 + v^ 2 - o; 2 r) . 



(3.6) 



i /-« 



("Mirror" Lifshitz) In this case the choice of the normalizable solution is unique 



and reads 



•u = ***■** i\-£- k .\>»* 



(3.7) 



We now use (3.1) to compute the correlators, up to normalization factors and contact 



terms. In the k — )■ limit the Schroedinger equation has the same form in all cases, with 



different coefficients in front of the 1/r 2 potential (2.10). The solutions are 



$(r) 



r 1+ lKi^i 



i\u\r 



(3.1 



where Kp(x) is the modified Bessel function of the second kind. In the cases we study we 
find that 

z = 2, A = 4, G 2 (lu,0) ~ co 2 log(inu), 
z = l, A = 3, G 2 (w,0) ~ |w| 3 , 
1 

2 



: = - A = |, G 2 (tu,0)~ \co\ 5 



(3.9) 



The behavior at non-zero momentum is more interesting. Let us now calculate the two-point 
correlators 



14 



z=2 (A = 4) 

G 2 {cu,k)~ 



Auj 2 



K' 



+ k* 



log (ikw) + ip ( — I + «G(Imo;)7rsech I — — 

2, Auj I \ Auj 



, (3.10) 



where ip(x) is the polygamma function and Q(x) is the Heaviside step function. The corre- 
lator has a branch cut along the positive imaginary axis. On top of it there are also poles at 
the positions 

7 K n 



There are two sources of poles in (3.10) - the polygamma function and the sech function. 



For n = 2m — 1 , m = 1, 2, . . . both functions contribute, for the other n only the ip has a 
pole giving a different residue at those poles. 

Notice that there is an infinite set of poles at any vicinity of u — 0. These poles all go to 
u = in the k — > limit. 



z = 1 (A = 3) 

G 2 (tu,k)~(k 2 -tu 2 ) 3/2 - 
This is the usual relativistic result, with a branch cut for values of u 2 > k 2 . 



(3.12) 



z = \ ( A = I) 



G 2 (u,k)~ k 5 ? 2 - 



r 



Z _ _y! 

4 2feft: 



im :5 — r (3 - 13) 

In this case the correlator is analytic in the complex frequency plane. There is a set of 
discrete poles localized on the real axis 

7 s 



and set of zeroes 



^',. 



w„ 



2n+-)/rife, n = 0,1,2,... 



2n--]K 2 A;, n = 0,l,2, ... 



(3.14) 



(3.15) 



When fc — > 0, the poles and the zeroes go to uj = 0. The poles of the correlator give a gapless 
spectrum of propagating modes, with a dispersion relation uj ~ yk. 

In order to show that these modes produce superluminal propagation we will follow 



the analysis in appendix C of ref. 25 . A wavefront is produced by a source of the form 



Q(t)e lut 5(x). In the linear approximation the response is approximately 



(0(x,t)) 



duo 

2^ 



^Rniu, k n ) 



-iu)t+ik n x 



uj — v + ie 



(3.16) 



where k n correspond to the modes (3.14) seen as poles in the complex momentum plane and 
R n are the residues. Specifically, 

,2 



k„ 



UJ 



k 2 (2n + I) ' 



n = 0,1,2, 



(3.17) 
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The phase velocity is defined as 

w c ( _ 7 



*>« = K = ^\ 2n+ 2)- (3 ' 18) 

The properly defined wavefront velocity is obtained by taking the large frequency limit of the 
phase velocity defined as above. Naively it seems that v p h would vanish, but notice that for 
any given large value of the frequency u> ^> £/c, there will always be an integer n w such that 
there is an infinite set of modes n > n^ with v p h > c. This implies that there is superluminal 
propagation, as we anticipated from the analysis of null geodesies. 

4 The NEC and Higher Derivative Gravity 

We have examined the relation between the null energy condition in the bulk and causality in 
the boundary theory for Lifshitz geometries and we have concluded that it should be satisfied 
in order to have a sensible holographic interpretation. A different kind of gravitational models 
that also produce Lifshitz geometries are higher derivative corrections of Einstein gravity, like 



the curvature squared corrections considered in ref. 26 27 . Higher derivative corrections 



of gravity and causality issues in holography have also been considered in the context of 



the AdS/CFT correspondence, mainly for Gauss-Bonnet gravity 28,29 , and also for quasi- 
topological gravity |30|. As we will see, our analysis is constrained to some particular class 
of models that in general do not include Gauss-Bonnet gravity, so we will not be able to 
impose further restrictions on this model. 

In principle our analysis only shows that the NEC is violated for particular solutions, 
determined by the values of A and the $ parameters. Notice that the cosmological constant 
A does not enter in the NEC, so that for the values of fa where a solution with z < 1 exists 
it is possible that the NEC is violated in more cases. For instance, Einstein gravity with a 
cosmological constant saturates the NEC, so if R 2 terms are introduced as small corrections, 
this implies that some mechanism should be at work to prevent general perturbations of 
the solutions from violating the NEC. A possibility could be to use boundary conditions 
such that problematic modes are projected out. This would require a thorough analysis of 
fluctuations, we will leave it as future work and take the simplest approach here. These 
arguments do not apply to Lifshitz solutions with z > 1, as long as fluctuations are small 
they should not spoil the NEC. Notice that there is no contradiction here, even if the R 2 
corrections are small, Lifshitz solutions have different asymptota to solutions to the Einstein 
equations, so they cannot be considered as small perturbations even if z is close to one. 



The authors of ref. 27 have made an extended analysis of various black hole solutions 
including Lifshitz black holes. The zero temperature Lifshitz metric is also a solution, and 
depending on the values of the parameters it is possible to find solutions with z > 1 or z < 1. 
Whenever a solution with z < 1 exists, this implies that the null energy condition can be 
violated for that choice of parameters. 

A general form of the action is 

S = J d D x^/g (R-2A + L 2 f3 x R 2 + L 2 f3 2 R a pR a ? + L 2 p 3 R a ^ s R a ^ s ) . (4.1) 
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The equations of motion read 



Rfj,u — yd^vR + Ag^ u — L <d^ u , 



(4.2) 



where Q^ stands for the variation of the higher derivative part of (4.1) 27 



+A/3 3 R lia R" — 2f3iRR ilv — 2(/3 2 + 2/3 s )R liau / 3 R a — 2/3 3 R lia i3 1 R\ 



pct/37 



(2 fa + /3 2 + 2fa)V lt V„R - -(4/5i + /3 2 )^V 2 J R - (fa + 4/3 3 )V 2 /^ 



(4.3) 



At the level of the equations of motion this tensor looks like some matter energy-momentum 



tensor that sources the Einstein equations. We can now look for Lifshitz solutions (1.5) of 



(4.2) which will partially fix the parameters A and /?, of the action (4.1). For a general 



number of dimensions D, there are Lifshitz solutions provided that 
1 



A 



1 + 2(fa - fa) + 2z+(l-2z + \z A ) (Afa + 2 fa + Afa 



+ (3z 2 -2z 3 )(f3 2 + 4f3 3 



(4.4) 



2(2z 2 + (D - 2)(2z + D - I)) fa + 2{z 2 +D- 2) fa + A(z 2 - (D - 2)z + I) fa = 1 • (4.5) 

The latter condition determines a plane in the (fa, fa, fa) parameter space for each value of 
z. Notice that all planes intersect at a single point 

fa = fa = -&/4 = 1/(2(Z> - A)(D - 3)) . (4.6) 

These values correspond to a particular case of Gauss-Bonnet gravity. Here and in the 



following we will assume that D > 4. If we solve for fa or /3 2 in (4.5) and plug the result in 



(44J), we find 

A: 



AD 



(2z + (D - 2)(2z + D - 1) - 4(L> - 3)(D - 4)z(z + D - 2)/3 3 



(4.7) 



There is a second branch of solutions with z = 1, where the cosmological constant is fixed 



in terms of the fa^s (4.4) but there is no constraint (4.5). Notice that the two branches will 



have the same value of parameters only when a z = 1 solution is allowed in the first branch, 
so the cosmological constant is the same. This happens at the points in the plane 

2D{D - I) fa + 2{D - I) fa - A(D - A) fa = 1 . (4.8) 

can be written as a homogeneous 



The condition for the existence of Lifshitz solutions (4.5 



equation for planes centered around the special point (4.6 



ai (z)(fa - fa) + a 2 (z)(fa - fa) + a 3 (z)(fa - fa) 







(4.9) 
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where 

D-l)), 



ai(z) = 


= 2(2z 2 + (L>-2)(2z + 


a 2 {z) -- 


= 2(z 2 + D-2), 


a 3 (z) = 


= 4(z 2 - (D- 2)z + l). 



(4.10) 

The set of planes defined by the values z G (—00, 00) span a volume in the (fa, fa, fa) 
space. Models with R 2 corrections to Einstein-Hilbert gravity where violations of the null 
energy condition is possible lie in the region — 00 < z < 1. Notice that the value of the 
cosmological constant is different for models with different values of z and the same values 
of the (3i coefficients. 

Given a value of z the theories that admit a Lifshitz solution lie on a plane that contains 
the origin of the Xi = (3i — $i space. The plane is determined by the equation a(z) • x = 0, 
where we use bold face for three-dimensional vectors. If we shift the value of the dynamical 
exponent z — > z + 5z, Lifshitz solutions will be allowed in the plane a(z + Sz) • x = 0. For 
an infinitesimal variation a(z + 8z) ~ a.(z) + a!(z)Sz the intersection between the two planes 
is given by the line 

a(z)-x = 0, a'0)-x = 0. (4.11) 

A general solution is given parametrically as 

Xi = T(a 2 (z)a' 3 (z)-a' 2 (z)a 3 {z)), 
X 2 = T(a3{z)a[ (z) - a' 3 (z)ai(z)), 
X 3 = r(oi(z)4(z)-oi(«)o 2 («)). (4.12) 

As we vary z, the lines span a surface, the envelope of the planes. The explicit form is 

Xi = 8r ((D - 2)z 2 + 2(D - 3)z - (D - 2) 2 ) , 

X 2 = -8t(A(D-2)z 2 + 2D(D-3)z-(D-2)(D(D-3)+A)), 

X 3 = 8t((D-2)z 2 + (D-2)(D-3)z-(D-2) 2 ). (4.13) 

Each plane is tangent to this surface along the line ( |4.11[ ). Furthermore, the equation a.(zi) ■ 
X(^ 2 , r) = is satisfied only if z\ = z 2 , so each plane intersects the surface X only once. We 
can solve z as a function of Xi as well, there are two possible solutions 

*=£r+£+2£ [(D - 2)(Ii -* ,)± 



±^{D - 2) 2 (xi - x 3 ) 2 - (2xi + x 2 + 2x 3 )((D -2){D- l)xi + (D - 2)x 2 + 2x 3 ) 

(4.14) 
The region with allowed Lifshitz solutions corresponds to values of Xi such that the argument 
of the square root is positive. The boundary of the allowed region is represented by the surface 

(D - 2) 2 ( Xl - S3) 2 - (2xi + x 2 + 2x 3 )((D - 2)(D - l)x x + (D - 2)x 2 + 2x 3 ) = 0. (4.15) 



This is just the implicit form of the envelope surface (4.13[), as one can easily check by 



introducing the explicit expressions in the equation (4.15) with Xi = X^. In order to avoid 



Lifshitz solutions with z < 1 we should impose z + > 1 and Z- > 1. This can be simplified if 
we impose first the condition z + z^ > 1, then discard negative values of z and finally discard 
the cases where 

z±> — >1. (4.16) 



We find that 



then 



Z + Z- 

((D-2)(D 
((D-2)(D 



[D - 2){D - l)xi + {D- 2)x 2 + 2x 3 



2xi + X2 + 2X; 



> 1 



2)xx 
2)xx 



(L>-3)x 2 >0 if 2xiH-x 2 + 2x 3 >0. 
(D-3)x 2 <0 if 2xiH-x 2 + 2x 3 < 0. 



and since z + and z_ have to be positive, we also have 

(D - 2)(D - l)x 1 + {D- 2)x 2 + 2x 3 > 
(D - 2){D - l)x! + {D- 2)x 2 + 2x 3 < 



if 
if 



2xi + x 2 + 2x 3 > 
2xi + x 2 + 2x 3 < 



(4.17) 



(4.18) 



(4.19) 



Notice that in this region the argument inside the square root in (4.14) is always smaller in 



absolute value than the first term squared, since (2xi + x 2 + 2x 3 )((.D — 2)(D — l)xi + (D 
2)x 2 + 2x 3 ) > 0. The regions where z± are positive are 



x 3 > xi if 
x 3 < xi if 



2xi + x 2 + 2x 3 > , 
2xi + x 2 + 2x 3 < . 



(4.20) 



As far as (4.16) is concerned, we can write z± 



B is the square root term in (4.14) and C 



= (A ± B)/C, where A = (D - 2)(x 3 - a*), 
2xi + x 2 + 2x 3 . From (4.20) we have the 



conditions that B > and A and C have the same sign. Let us first consider the case C > 0, 
A > 0. Clearly z + > z_, now imagine that we are in the situation we want to discard, when 
z + > 1 > Z-. We can do the following manipulations 



Z- 



A+B 

£b 

c 



>1, 

< 1, 



A + B>C 
A-B <C 



B>0. 



(4.21) 



We know that the last expression is true, so the only allowed possibility is z+ — Z- — 1. The 
case C < 0, A < works in a similar way, if we consider the case Z- > 1 > z+, 



Z- 



A+B 

£b 

c 



<1, 

>1, 



A + E^C 
A-B <C 



B > 



(4.22) 



which again is always true. Therefore the NEC can be violated in the full region of param- 
eter space where Lifshitz solutions with z ^ 1 exist. Our results for the case D = 5 are 
summarized in Fig. [3j 

5 Conclusions and Open Questions 

In this paper we have investigated the role of the null energy condition in a gravitational 
model with respect to causality in holography. We have shown that the NEC has a simple 
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ih »■» 




Figure 3: Solutions with Lifshitz scaling are allowed in the colored regions of (/3i, /?2, Pa) space (left). 
Solutions with z < 1 exist in the full region for determined values of the cosmological constant, 
so violations of the NEC are possible in the full region. The thick dot in the origin corresponds 
to Einstein gravity. The right plot is a slice on the Pi = plane, where the dot corresponds to 
Einstein gravity, we observe that it lies at the boundary of the Lifshitz region, where two different 
z = 1 solutions exist. 



physical manifestation in terms of the local speed of light in the bulk. If the local speed 
of light on a radial slice increases when an observer moves towards the boundary, then and 
only then the NEC is satisfied. If we consider the Lifshitz metric with the critical exponent 



z (1.5), the NEC merely requires z > 1. The behavior of the speed of light when the NEC is 



violated implies that shock wave sources can bounce back to the boundary and reach it at 
points that are outside the light cone of signals propagating along the boundary. The front of 
radiation produced by the shock wave will also move at superluminal speeds on the average. 
Based on the results of ref. 
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we have argued that this would lead to superluminal modes, 
so we have continued with the investigation of the spectra of scalar field perturbations in 
Lifshitz backgrounds for different values of the critical exponent. We have shown explicitly 
that for a Lifshitz background with z — 1/2 there is indeed a set of discrete modes whose 
dispersion relation is incompatible with causality, while for examples with z > 1 such modes 
are absent. We have related the presence of such modes to the properties of the equation of 
motion for backgrounds with z < 1, by showing that the equation written in the Schrodinger 
form has a confining potential at non-zero momentum. 

Summarizing, we have presented strong arguments and some evidence that geometries 
produced by matter that violates the NEC will produce superluminal propagation in the 
dual theory. This is directly related to the observation in higher derivative theories that 
fluctuations in AdS with negative energy fluxes will spoil causality in the CFT 28 ,j It 



3 A recent work also connects negative energy fluxes to the presence of ghosts in the CFT 



31 
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would be interesting to confirm this statement in more general cases. There is also still an 
open question of whether violations of the NEC in a localized region of the bulk, like those 
produced by quantum effects, are allowed in holographic systems. It is also interesting to 
observe that the NEC seems to be an important condition for several aspects of holography, 
in principle not directly related to each other in the field theory. In addition to causality 



constraints, it allows the formulation of holographic c-theorems 22 and it was originally 



proposed in the context of the entropy bound 32 . A possible application of these results 
is to constrain extensions of holographic models with new forms of matter or curvature 
corrections to Einstein gravity. We have done this for holographic models that involve 
curvature squared corrections to Einstein gravity. Our criterion was to discard the regions 
of the parameter space that allow Lifshitz solutions with z < 1 and hence violations of the 
NEC by small perturbations around the solutions to Einstein equations. More generally, one 
would constrain possible holographic constructions by separating Einstein equations in two 
parts - one being the Einstein tensor and the other an effective energy-momentum tensor 
that includes the higher curvature corrections and any additional matter and imposing the 
null energy condition on the last. 

We have worked mainly in the bulk theory, let us now draw our attention to the macro- 
scopical properties of the boundary theory. We can consider introducing a finite temperature 



in the Lifshitz theory, or in the dual description a black hole 33 . Given the scale invariance 



of the theory, the equation of state for a perfect fluid should have the form 

z(T tt )-d(T xx )=0. (5.1) 

Assuming the metric of the boundary theory is flat let us consider the so-called dominant 
energy condition (DEC) for the boundary theorjF] The DEC is responsible for the causal 
flow of the energy-momentum. In covariant formulation is states 

-T°r>0, (5.2) 

for any light-like or timelike vector ^ such that £° > 0. For the theory with spatial isotropy 
the DEC reads 

(T tt ) > 0, (T tt ) - \(T XX )\ = (T tt ) (l - £) > 0, (5.3) 



where we used (5.1). The boundary DEC is satisfied provided z < d. This condition, 
combined together with the constraint from causality we have discussed in Sec. [2] gives the 
domain of critical exponents 

\<z<d. (5.4) 

Notice that the upper bound has been derived from the boundary DEC, while the lower 
bound has been derived from the bulk NEC. It is not clear how seriously the upper bound on 



the critical exponent should be considered. For instance, the supergravity construction of 34 



allows for a Lifshitz background with z ~ 39; even real condensed matter systems 35 with 
z > d are known. The most likely explanation, at least for the condensed matter systems, is 
that the DEC cannot be applied to the energy-momentum tensor we are considering. 



7 Boundary NEC is automatically satisfied for positive energy density. 
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